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(Dated: July 5, 2001) 
We present an updated measurement of time-dependent CP- violating asymmetries in neutral 
B decays with the BABAR detector at the PEP-II asymmetric B Factory at SLAC. This result 
uses an additional sample of T(4S') decays collected in 2001, bringing the data available to 32 
million BB pairs. We select events in which one neutral B meson is fully reconstructed in a 
final state containing charmonium and the flavor of the other neutral B meson is determined 
from its decay products. The amplitude of the CP-violating asymmetry, which in the Standard 
Model is proportional to sin2/3, is derived from the decay time distributions in such events. The 
result sin2/3 — 0.59 ± 0.14 (stat) ± 0.05 (syst) establishes CP violation in the S" meson system. 
We also determine |A| = 0.93 ± 0.09 (stat) ± 0.03 (syst), consistent with no direct CP violation. 

PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 



CP violation has been a central concern of particle 
physics since its discovery in 1964 in the decays of K^ 
mesons |1|. To date, this phenomenon has not been ob- 
served in any other system. An elegant explanation of 
this effect was proposed by Kobayashi and Maskawa, as a 
CP-violating phase in the three-generation CKM quark- 
mixing matrix |3]. In this picture, measurements of CP- 
violating asymmetries in the time distributions of B^ de- 
cays to charmonium final states provide a direct test of 
the Standard Model of electroweak interactions B, free 
of corrections from strong interactions that obscure the 
theoretical interpretation of the observed CP violation in 
K^ decays. 

Measurements of the CP-violating asymmetry param- 
eter sin2/3 have recently been reported by the BABMl Q] 
and Belle collaborations, from data taken in 1999- 

2000 at the PEP-II and KEKB asymmetric-energy e+e~ 
colliders respectively, with better precision than previous 
experiments n . In this Letter we report a new measure- 
ment of sin2/3, enhanced by 9 million BB pairs collected 
in 2001, additional decay modes, and improvements in 
data reconstruction and analysis. The BABAR detector 
and the experimental method are described in Refs. |j, Q, 
so the discussion here is limited to items and issues per- 
tinent to the current analysis. 

The complete data set (32 million BB pairs) has 
been used to fully reconstruct a sample Bcp of neutral 
B mesons decaying to the J/4>K°, tI>{2S)K°, J/iI>KI, 
XciK°, and J/iPK*°{K*° -^ K°tt°) final states. The 
last two modes have been added since Ref. M. There 
are several other significant changes in the analysis. Im- 
provements in track and K'^ reconstruction efficiency in 

2001 data produce an approximately 30% increase in the 
yields for a given luminosity. Better alignment of the 
tracking systems in 2001 data and improvements in the 
tagging vertex reconstruction algorithm increase the sen- 
sitivity of the measurement by an additional 10%. Opti- 
mization of the J/^IjK^ selection increases the purity of 
this sample. The final Bcp sample contains about 640 
signal events and, with all the improvements, the statis- 



tical power of the analysis is almost doubled with respect 
to that of Ref. g. 

We examine each of the events in the Bcp sample for 
evidence that the other neutral B meson decayed as a 
B^ or a B^ (fiavor tag). The decay-time distributions 
for events with a i?° or a B^ tag can be expressed in 
terms of a complex parameter A that depends on both 
j^oj^o Ynbimg and on the amplitudes describing B*^ and 
B^ decay to a common final state / m . The distribution 
f+(f_) of the decay rate when the tagging meson is a 
B°{B^) is given by 
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-lAtl/rgo 



2tbo(H-|A|2) 
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, (1) 



where At = tcp — ^tag is the time between the two B 
decays, r^o is the i?° lifetime and ArriBO is the mass dif- 
ference determined from B^B^ mixing |^ . The first oscil- 
latory term in Eq. |l| is due to interference between direct 
decay and decay after mixing. A difference between the 
B'^ and S" At distributions or a At asymmetry for either 
fiavor tag is evidence for CP violation. 

In the Standard Model A = rjfe"'^'^^ for charmonium- 
containing b -^ ccs decays, rjf is the CP eigenvalue of 
the state / and /3 = arg [ -KdKb/Kd^tb ] i^ ^^ angle of 
the Unitarity Triangle of the three-generation CKM ma- 
trix 0. Thus, the time-dependent CP-violating asym- 
metry is 



AcpiAt)^'-^^'^-'-^^'^ 



f+(Ai) + f-(Ai) 
= —rjf sm2f3 sin (Amso At) , 



(2) 



where tjf = ~1 for J/tpK^, ij{2S)K° and XciK° and +1 
for J/ipKl. Due to the presence of even (L=0, 2) and odd 
(L=l) orbital angular momenta in the J/4>K*'^{K*'^ — > 
A'gTr") system, there can be CP-even and CP-odd contri- 
butions to the decay rate. When the angular information 
in the decay is ignored, the measured CP asymmetry in 



320 



^240 h 

r^i6o 



C 80 



I I J/\|(K^(3I*7l") 

E3 J/\(fKg(7cV) 

□ J/\|(K*°(K°Jl") 

^ \(((2S)K^(7lV) 



n 



a) 



^ 

5.2 




5.22 



5.24 



5.26 



5.28 



5.3 



mp^ (GeV/c') 




4 £)ato 
I I J/\if Kl signal MC 
^^ Background 



't » * *« ;:*^FP* 



b) 



60 80 

AE (MeV) 

FIG. 1; a) Distribution of tties for Bcp candidates having 
a Ks in the final state; b) distribution of IS.E for J/i/'i^'^ 
candidates. 



J/'ipK*'^ is reduced by a dilution factor D± = 1 — 2i?j^, 
where R± is the fraction of the L=l component. We have 
measured R±_ = (16 ±3.5)% Q which, after acceptance 
corrections, leads to an effective rjf = 0.65 ± 0.07 for the 
J/il:K*° mode. 

The hadronic event selection, lepton and charged kaon 
identification, and J/f/; and i}}{2S) reconstruction rele- 
vant to this analysis have been described in Ref. M, as 
have the selection criteria for the channels J/ijj K^ {K^ -^ 
7r+7r-,7r°7r°), ^|^{2S)K° {K° -^ Tr+Tr"), and J/^pK1. In 
the J/tp Kl selection, the transverse missing momentum 
requirement has been reoptimized for the Acp study. 

For the decay i3° -^ XciK'^, the mode Xci ~* J/''!'! 
is reconstructed with mass-constrained J/ip candidates 
selected as in other charmonium channels Q. Photons 
must have an energy greater than 150 MeV and must not 
be associated with any reconstructed tt*'. The resulting 
J/tpj mass is required to be within 35 MeV/c^ of the Xci 
mass p. 

For the decay B° -^ J/%IjK*^, the K*^ -^ K°7r° candi- 
date is formed by combining a 7r° — > 77 candidate satisfy- 
ing 106 < m-yj < 153 MeV/ c^ with a K'^ candidate. The 
cosine of the angle between the K'^ momentum vector in 
the K*'^ rest frame and the K*'^ momentum defined in 
the B rest frame is required to be less than 0.95. We 
require 796 < m^^o^o < 996MeV/c2. 

Bcp candidates are selected by requiring that the dif- 
ference AE between their energy and the beam energy 
in the center-of-mass frame be less than 3cr from zero. 
For modes involving K'^, the beam-energy substituted 



mass mES = \/(-^bGam)^ ~ (P™)^ must be greater than 
5.2GeV/c^. The resolution for AE is about 10 MeV, ex- 
cept for the if° -^ ttOtt" mode (33 MeV), the J/?/'if*° 
(20 MeV) and the J/i/^Kl modes (3.5 MeV after B mass 
constraint). For the purpose of determining numbers of 
events and purities, a signal region mes > 5.27 GeV/c^ is 
used for all modes except J/ipK^ and J/tpK*'^. 

Figure m shows the resulting toes distributions for Bcp 
candidates containing a K'^ and the AE distribution for 
the candidates containing a K^. The Bcp sample con- 
tains 1230 events in the signal region (before tag and 
vertex requirements), with an estimated background of 
200 events, predominantly in the J/ipK^ channel. For 
that channel, the composition, effective 77/, and AE dis- 
tributions of the individual background sources are taken 
either from a Monte Carlo simulation (for B decays to 
J/tp) or from the rrii+g- sidebands in data. 

A measurement of Aqp requires a determination of the 
experimental Af resolution and the fraction of events in 
which the tag assignment is incorrect. A mistag fraction 
w reduces the observed asymmetry by a factor (1 — iw). 
A sample of self-tagging B decays -Bflav used in the deter- 
mination of the mistag fractions and At resolution func- 
tions consists of the channels D(*-'^ft,+ (/i+ = 7r+, p+, Oj^) 
and J/^K*° {K*° -^ K+tt-) 0. A control sample of 
charged B mesons decaying to the final states J/tpK^*^^, 
Tp{2S)K+, XciK+ and i5(*)"7r+ is used for validation 
studies. 

For flavor tagging, we exploit information from the 
other B decay in the event. Each event is assigned to 
one of four hierarchical, mutually exclusive tagging cat- 
egories or excluded from further analysis. The Lepton 
and Kaon categories contain events with high momen- 
tum leptons from semileptonic B decays or with kaons 
whose charge is correlated with the flavor of the decay- 
ing b quark {e.g. a positive lepton or kaon yields a B'^ 
tag). The NTl and NT2 categories are based on a neural 
network algorithm whose tagging power arises primarily 
from soft pions from D*^ decays and from recovering 
unidentified isolated primary leptons ^. 

The numbers of tagged events and the signal purities, 
determined from fits to the toes (all K'^ modes except 
K*^) or AE {K^ mode) distributions in data or from 
Monte Carlo simulation {K*'^ mode), are shown in Ta- 
ble |. The efficiencies and mistag fractions for the four 
tagging categories are measured from data and summa- 
rized in Table |l|. 

The time interval At between the two B decays is then 
determined from the Az = zcp — -Ztag measurement, in- 
cluding an event-by-event correction for the direction of 
the B with respect to the z direction in the T{4S) frame. 
zcp is determined from the charged tracks that constitute 
the Bcp candidate. The tagging vertex is determined by 
fitting the tracks not belonging to the Bcp (or Bflav) 
candidate to a common vertex. The method employed is 
identical to our previous analysis except for the addition 



TABLE I; Number of tagged events, signal purity and result 
of fitting for CP asymmetries in the full CP sample and in 
various subsamples, as well as in the -Bflav and charged B 
control samples. Errors are statistical only. 



A^tag Purity (%) 



sin2/3 



Sample 

J/ip K'i,'ipi2S)K'i,XciK'i 480 96 0.56 ± 0.15 

J/ip Kl {rif ^ +1) 273 51 0.70 ± 0.34 

J/tPK*°,K*° -^ K°s-k° 50 74 0.82 ± 1.00 

Full CP sample 



803 



80 



0.59 ± 0.14 



J/^K%, H2S)K)i, XciK^s only fa" 



JI^K"s {K^s^ 


TT+TT") 


316 


98 


0.45 ± 0.18 


J/^K° {K^s^ 


TV TT ) 


64 


94 


0.70 ± 0.50 


^{2S)Kt [K% 


~*7V + TT~) 


67 


98 


0.47 ± 0.42 


XciKl 




33 


97 


2.59 ± []:i? 


Lepton tags 




74 


100 


0.54 ± 0.29 


Kaon tags 




271 


98 


0.59 ± 0.20 


NTl tags 




46 


97 


0.67 ± 0.45 


NT2 tags 




89 


95 


0.10 ± 0.74 



B" tags 234 
B° tags 246 


98 
97 


0.50 ± 0.22 
0.61 ± 0.22 


Bflav non-CP sample 7591 


86 


0.02 ± 0.04 


Charged B non-CP sample 6814 


86 


0.03 ± 0.04 



TABLE II: Average mistag fractions Wi and mistag differ- 
ences Awi — Wi{B ) — Wi{B°) extracted for each tagging 
category i from the maximum-likelihood fit to the time dis- 
tribution for the fully-reconstructed B" sample (Bflav +Bcp). 
The figure of merit for tagging is the effective tagging effi- 
ciency Qi — £i{l — 2'Wif' , where Si is the fraction of events 
with a reconstructed tag vertex that are assigned to the i"* 
category. Uncertainties are statistical only. The statistical 
error on sin2/9 is proportional to l/\/Q, where Q — ^ Qi. 



Category 



Lepton 
Kaon 
NTl 
NT2 



:(%) 



w 



10.9 ± 0.3 

35.8 ± 0.5 

7.8 ± 0.3 

13.8 ± 0.3 



(%) 



8.9 ±1.3 
17.6 ±1.0 

22.0 ±2.1 

35.1 ±1.9 



Aw (%) 
0.9 ± 2.2 



Q(%) 

7.4 ± 0.5 



-1.9 ± 1.5 15.0 ± 0.9 

5.6 ± 3.2 2.5 ± 0.4 

-5.9 ± 2.7 1.2 ± 0.3 



All 



;.4 ± 0.7 



26.1 ± 1.2 



of a constraint from knowledge of the beam spot loca- 
tion and beam direction, which increases its efficiency 
to 97%. This is incorporated through the addition of a 
pseudotrack to the tagging vertex, computed from the 
Bqp (Bflav) vertex and three-momentum, the beam spot 
(with a vertical size of 10 /im) and the T{AS) momentum. 
For 99% of the reconstructed vertices the r.m.s. Az res- 
olution is 180 /im. An accepted candidate must have a 
converged fit for the Bcp and i^tag vertices, an error of 
less than 400 //m on Az, and a measured |Ai| < 20 ps. 
After tag and vertexing requirements about 640 signal 
events remain. 

The sin2/3 measurement is made with a simultaneous 
unbinned maximum likelihood fit to the At distributions 
of the Bcp and Bflav tagged samples. The At distribu- 
tion of the former is given by Eq. |l|, with |A| = 1. The 



Bflav sample evolves according to the known rate for fla- 
vor oscillations in neutral B mesons H. The amplitudes 
for Bcp asymmetries and for Bflav flavor oscillations are 
reduced by the same factor (1 — 2w) due to wrong tags. 
Both distributions are convolved with a common At res- 
olution function and corrected for backgrounds. Events 
are assigned signal and background probabilities based 
on the mES (all modes except J/ipK^) or AE {J/ipK^) 
distributions. 

The representation of the At resolution function is the 
same as in B with small changes: all offsets are modeled 
to be proportional to crAt, which is correlated with the 
weight that the daughters of long-lived charm particles 
have in the tag vertex reconstruction. Separate resolution 
functions have been used for the data collected in 1999- 
2000 and 2001, due to the significant improvement in the 
silicon vertex tracker (SVT) alignment. The scale factor 
for the tail component is fixed to the Monte Carlo value 
since it is strongly correlated with the other resolution 
function parameters. 

A total of 45 parameters are varied in the likelihood 
fit, including sin2/3 (1), the average mistag fraction w 
and the difference Aw between B° and B" mistags for 
each tagging category (8), parameters for the signal At 
resolution (16), and parameters for background time de- 
pendence (9), At resolution (3) and mistag fractions (8). 
The determination of the mistag fractions and signal At 
resolution function is dominated by the large Bgav sam- 
ple. Background parameters are governed by events with 
Toes < 5.27GeV/c^. As a result, the largest correlation 
between sin2/3 and any linear combination of the other 
free parameters is only 0.13. We fix t^o = 1.548 ps and 
Atobo = 0.472 ?ips-i §. 

Figure shows the At distributions and Acp as a func- 
tion of At overlaid with the likelihood fit result for the 
rjf = —\ and r// = +1 samples. The probability of ob- 
taining a lower likelihood, evaluated using a Monte Carlo 
technique, is 27%. The simultaneous fit to all CP decay 
modes and fiavor decay modes yields 

sin2/3 = 0.59 ± 0.14 (stat) ± 0.05 (syst). 

Repeating the fit with all parameters fixed to their deter- 
mined values except sin2/3, we find a total contribution 
of ±0.02 to the error on sin2/3 due to the combined sta- 
tistical uncertainties in mistag fractions. At resolution 
and background parameters. The dominant sources of 
systematic error are the parameterization of the At reso- 
lution function (0.03), due in part to residual uncertain- 
ties in SVT alignment, possible differences in the mistag 
fractions between the Bcp and Bflav samples (0.03), and 
uncertainties in the level, composition, and CP asymme- 
try of the background in the selected CP events (0.02). 
The systematic errors from uncertainties in Attibo and 
Tga and from the parameterization of the background in 
the Bflav sample are small; an increase of 0.02 Tips^^ in 
the value for Ato^o decreases sin2/3 by 0.015. 
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FIG. 2: Number of??/ = -1 candidates (J/V'A's, ■(/'(2S')iC2, 
and Xci^s) in the signal region a) with a B" tag A^^o 
and b) with a -B" tag Ngo, and c) the asymmetry (Ngo — 
^B'^)/i^BO +Ngo), as functions of At. The solid curves rep- 
resent the result of the combined fit to all selected CP events; 
the shaded regions represent the background contributions. 
Figures d)-f) contain the corresponding information for the 
r)f = +1 mode [J/tpKl). The likelihood is normalized to the 
total number of B° and B" tags. The value of sin2/3 is inde- 
pendent of the individual normalizations and therefore of the 
difference between the number of S" and B'^ tags. 



The large sample of reconstructed events allows a num- 
ber of consistency checks, including separation of the 
data by decay mode, tagging category and Stag flavor. 
The results of fits to these subsamples are shown in Ta- 
ble |. The consistency between various modes is satisfac- 
tory, the probability of finding a worse agreement being 
8%. The observed asymmetry in the number of B^ (160) 
and B^ (113) tags in the J/ipK^ sample has no impact 
on the sin2/3 measurement. Table | also shows results of 
fits to the samples of non-CP decay modes, where no sta- 




tistically significant asymmetry is found. Performing the 
current analysis on the previously published data sam- 
ple and decay modes yields a value of sin2/3 =0.32±0.18, 
consistent with the published value [Q . For only these de- 
cay modes, the year 2001 data yield sin2/3 =0.83±0.23, 
consistent with the 1999-2000 results at the 1.8a level. 

If |A| is allowed to float in the fit to the 77/ = —1 sample, 
which has high purity and requires minimal assumptions 
on the effect of backgrounds, the value obtained is |A| = 
0.93 ± 0.09 ((stat)) ± 0.03 ((syst)). The sources of the 
systematic error in this measurement are the same as in 
the sin2/3 analysis. The coefficient of the sin(AmBoAi) 
term in Eq. ^ is measured to be 0.56 ± 0.15 (stat). 

The measurement of sin2/3 = 0.59 ± 0.14 (stat) ± 
0.05 (syst) reported here establishes CP violation in the 
B^ meson system at the 4.1<t level. This significance 
is computed from the sum in quadrature of the statis- 
tical and additive systematic errors. The probability of 
obtaining this value or higher in the absence of CP vi- 
olation is less than 3 x 10^^. The corresponding prob- 
ability for the rjf = —1 modes alone is 2 x 10^**. This 
direct measurement is in agreement with the range im- 
plied by measurements and theoretical estimates of the 
magnitudes of CKM matrix elements [|lj| . 
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